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Abstract. A differential-algebraic approach to studying the Lax type integrability of the generalized 
Riemann type hydrodynamic hierarchy is revisited, its new Lax type representation is constructed 
in exact form. The related bi-Hamiltonian integrability and compatible Poissonian structures of the 
generalized Riemann type hierarchy are also discussed by means of the gradient-holonomic and geometric 
methods. 



1. Introduction 

Recently new mathematical approaches, based on differential-algebraic and differential geometric 
methods and techniques, were applied in works [TJ I15( 1 1 7 j for studying the Lax type integrability of 
nonlinear differential equations of Korteweg-de Vries and Riemann type. In particular, a great deal of 
analytical studies [2J [3J [51 [TJ [5] were devoted to finding the corresponding Lax-type representations of 
the infinite Riemann type hydrodynamical hierarchy 

(1.1) D?u = 0, D t :=d/dt + ud/dx, 

where N G Z+, (x,t) T G R 2 and u G C°°(R/27rZ; R). It was found that the related dynamical system 

(1.2) D t ui=U2, -., D t Uj = u j+ i, ... ,D t u N = 0, 

defined on a 27r-periodic infinite-dimensional smooth functional manifold Mn C C°°(R/27rZ; R ), pos- 
sesses [3j [17] for an arbitrary integer N G Z + a suitable Lax type representation 

(1.3) D x f = l N [u;X]f, D t f = q N (\)f 

with A £ C being a complex spectral parameter and / G L !X ,(R;C' V ) and matrices In[u; A], <7jv(A) G 
EndC 2 . Here, by definition, u% := u G C°°(R 2 ;R) and the differentiations 

(1.4) D t := d/dt + u x D x , D x := d/dx 
satisfy on the manifold M/v the following commutation relationship: 

(1.5) [D x ,D t ) = (D x u 1 )D x , 
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In particular, for the cases N — 2,3 and N = 4 the following exact matrix expressions 
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polynomial in A € C, were presented in exact form. A similar Lax type representation, as it follows from 
the statements of the work [T7], also holds for the general case N G Z+. These results are mainly based 
on a new differential-algebraic approach, devised recently in work |17j and contain complicated enough 
[3] functional expressions r$ , j — 1, N — 2, which satisfy the recurrent differential- functional equations 



(1.7) 



f Ar«+rW^W = l, 



u t' N 



on the functional manifold Mat. Recently Popowicz [3] has found a new Lax representation for the 
generalized Riemann equation (|1.2[) at arbitrary A 6 Z + . In the present article, being strongly based on 
the methods devised in works [T71 [3] , we will revisit the differential-algebraic approach to the Riemann 
type hydrodynamical hierarchy (jl.l[) and construct its new Lax type representation 
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in a very simple and useful for applications form for any N € Z+ , which is scale equivalent to that found 
by Popowicz [4]. Moreover, we prove that this Riemann type hydrodynamical hierarchy generates the 
bi-Hamiltonian flows on the manifold Mjv and, finally, we analyze for the cases N — 2,3 and N = 4 
the corresponding compatible [23, 25, 27. Poissonian structures, following within the gradient-holonomic 
scheme from these new Lax type representations. A mathematical nature of the Lax type representations 
(|1.3j) . (|1.6[) presents from the differential- algebraic point of view a very interesting question, an answer 
to which may be useful for the integrability theory, remains open and needs additional investigations. 
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It is also worth to mention that the methods devised in this work can be successfully applied to other 
interesting for applications 16, 20, 2l][T9l[22] nonlinear dynamical systems, such as Burgers, Korteweg-de 
Vries and Ostrovsky-Vakhnenko equations 

(1.9) D t u = Dlu, D t u = D\u, D x D t u = -u, 

on the 27r-periodic functional manifold M\ C C^ 00 ) (R/27rZ; R), new infinite hierarchies [BJ|!|] of Riemann 
type hydrodynamic systems 

(1.10) D?- 1 u = D*z, D t z = Q, 
and 

(1.11) D^- l u = zl, D t z = 

on a smooth 27r-periodic functional manifold Mjv C C^ 00 ' (R/27rZ; R w ) for N £ N and many others. 

2. DIFFERENTIAL-ALGEBRAIC APPROACH REVISITING 

We will consider the ring K, := R{{a;,i}}, (x,ty £ R 2 , of the convergent germs of real-valued smooth 
functions from C°° (R 2 ; R) and construct [TU1 EJ [HI EH 03] the associated differential polynomial ring 
K{u} :— /C[0u] with respect to a functional variable u± := u £ JC, where O denotes the standard monoid 
of the all commuting differentiations d/dx and d/dt. The ideal I{u} C K.{u} is called differential if the 
condition I{u} = <dl{u} holds. 

In the differential ring IC{u}, interpreted as an invariant differential ideal in K., there arc naturally 
defined two differentiations 

(2.1) D t , D x : K{u} -> K{u}, 

satisfying the Lie-commutator relationship ()1.5j) . For a general function u £ K, there exists the only 
representation of (|1.5|) in the ideal JC{u} of form (|1.4|) . Nonetheless, if the additional constraint (jl.ljl 
holds, its linear finite dimensional matrix representation in the corresponding functional vector space 
JC{u} N for N £ Z+, related with some finitely generated invariant differential ideal I{u} C IC{u}, 
may exists being thereby equivalent to the corresponding Lax type representation for the Riemann type 
dynamical system (jl.2[) . To make this scheme analytically feasible, we consider in detail the cases TV = 1, 4 
and construct the corresponding linear finite dimensional matrix representations of the Lie-commutator 
relationship (|1.5|) . polynomially depending on an arbitrary spectral parameter A £ C. 

2.1. The case N = 1. Aiming to find the corresponding representation vector space for the Lie-algebraic 
relationship (|1.5j) at N = 1, we need firstly to construct |17j a so called invariant generating Riemann 
differential ideal R{u} C IC{u} as 

(2.2) R{u} := { E E ^ {n+l) fk 3+l) DlD n x u £ K{u} : /^ +1 > £ K, J, n £ Z + }, 



where A £ R\{0} is an arbitrary parameter. The differential ideal (I2.2[) is, evidently, invariant and 
characterized by the following lemma. 

Lemma 2.1. The kernel KerD t C R{u} of the differentiation D t : K.{u] — > K.{u] is generated by 
elements f^> £ R{{x,t}},j £ Z + , satisfying the linear differential-functional relationships 

(2.3) = -A/ W , 
where, by definition, 

(2.4) / (3+1) (A) = ^ /^ +1) A" n 
for A G R\{0} and j £ Z + . 



YAREMA A. PRYKARPATSKY 1 , OREST D. ARTEMOVYCH 2 , MAXIM V. PAVLOV 3 AND ANATOLIY K. PRYKARPATSKY 4 

Taking now into account the invariant reduction of the differential ideal (|2.2p subject to the condition 
D t u = to the ideal 

(2.5) RW{u] := { £ X ~ n fL 1)D > e : A« = 0, € AC, n € Z + } 

one can simultaneously reduce the hierarchy of relationships (12 .3[) to the simple expression 



(2.6) D 4 /«(A)=0, 

where Z' 1 - 1 £ /Ci{u} := /C{u}}|r> t ti=o- Now, to formulate the Lax type integrability condition for the case 
N = 1, it is necessary to construct the corresponding _D t -invariant Lax differential ideal 

(2.7) L«M := { 5 i/ (1) (A) G /CrH : A/ (1) (A) = 0, ffl 6 /c} 

and to check its invariance subject to the differentiation D x : ]Ci{u} — > K,\{u\. Namely, the following 
lemma holds. 

Lemma 2.2. The Lax differential ideal \2.T§ is D x -invariant, if the linear equality 

(2.8) A/ (1) (A) = (Xu x + tiu- 2 u xx )fW(\) 

holds for arbitrary parameters A,/i€l, where the subscript "x" means the usual D x - differentiation. 



Proof. The condition (|2.8|) makes the Lax ideal (|2.7|) to be also D^-invariant if the following condition 
on the linear representation of the Z^-differentiation 

(2.9) D w fW(A) = h[u;\]fW(\) 
holds: 

(2.10) D t h[u; X] + h[u; X]D X U! =0 

for any parameter A 6 R. Differential-functional equation (|2.10j) upon substitution 

(2.11) h[u;X]:=D xai [u;X\ 
reduces easily to the equation 

(2.12) D t a 1 [u;X] = X 

for any parameter A € R. Taking into account the condition D t u\ = one can easily hnd that 

(2.13) ai[u; A] = Ax + n/D x u\, 

where \i £ R is arbitrary. Having now substituted expression (|2.13[) into (|2 . 1 1 1) and taking into account 
that D t x = Mi, the result (12.80 follows. □ 



Thereby, having naturally extended the differential ring K, to the ring C{{x,i}}, one can formulate 
the following proposition. 

Proposition 2.3. The Lax type representation for the Riemann hydrodynamical equation 

(2.14) D t u = 

is given by a set of the linear compatible equations 

(2.15) A/ (1) (A) = 0, D x fW(X) = (Xu x +im- 2 u xx )fW(X) 
/or/W(A) eLoo(R 2 ;C) and for any A, [i € C. 
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2.2. The case N = 2. For the case N = 2 the respectively reduced invariant Riemann differential ideal 
(|2.2|) is given by the set 

R^{u} : \-^f^DiD«u:Dlu = 0, 

j=0,l nGZ + 

(2.16) g /C, ngZ+, i = M} 

and characterized by the kernel KerD t C i?( 2 ){u} of the D t -differentiation, which is generated by the 
following linear differential relationships: 

(2.17) A/ (1) (A) = 0, A/ (2) (A) = -A/«(A), 
where /«(A) g /C 2 {w} := K{u}\ 

_D 2 u=0: J = Oj 1; an d A g R is an arbitrary parameter. 
The condition (|2~T7f can be rewritten in a compact vector form as 



(2-18) £> t /(A) = <Z 2 (A)/(A), <? 2 (A) 





-A 



(2.21) l 2 [u;X] = 



where /(A) := (/«(A), / (2) (A)) T £ M"} 2 - 

Now we can construct the invariant Lax differential ideal 

2 

(2.19) L^{u} : = {»/ °(A) £ /C 2 { M } : A/ (1) (A) = 0, 

i=i 

A/ (2) (A) = -A/W(A),. 9j g/C, j = 172} 
and to check its invariance subject to a linear representation of the D^-differentiation in the form: 

(2.20) D x f(X)=l 2 [u;X]f(X) 

for some matrix l 2 [u; A] g End R 2 . 
The following lemma holds. 

Lemma 2.4. The Lax differential ideal i2.19\) is D x -invariant if the matrix 

Xui yX U2,x 

— 2A — Awi j: 

Proof. The D^- invariant condition for the Lax ideal (I2.19[) forces the matrix l 2 [u] A] g End R 2 to satisfy 
the differential-functional relationship 

(2.22) D t l 2 +l 2 D x ui = [q 2 ,l 2 ], 

where we have put by definition l 2 ■= l 2 [u; A], q 2 := 92(A) g End R 2 . Making the substitution 

(2.23) h = D x a 2 

for some matrix a 2 '■= a 2 [u; A] g End R 2 , we can easily reduce equation (|2.22[) to the equivalent one in 
the form 

(2.24) D t a 2 = [q 2 ,a 2 ]. 
To solve equation (|2.24l) it is useful to take into account that 

(2.25) D\a 2 = k 2 q 2 , D*a 2 = 

for some constant k 2 g R. Taking into account the relationship (I2.25|) one easily obtains the exact matrix 
representation 

, / Aui u 2 

(2 - 26) ° 2 = { -2X 2 x -X Ul 

entailing the result (|2~2l) □ 

Thus, having as above extended the differential ring K, to the ring C{{ir,i}}, one can formulate the 
next proposition. 
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Proposition 2.5. The Lax representation for the Riemann type hydrodynamical system 

(2.27) D tUl = u 2 , D t u 2 = 
is given by a set of the linear compatible equations 

(2.28) A/(A)=« 2 (A):=(_° A ° ) /(A), D x f(X) = ( ^ ) /(A) 
/or /(A) G Loo(IR 2 ;C 2 ) and arbitrary complex parameter A G C. 

2.3. The case N = 3. Similarly to the above, for the case N = 3 the respectively reduced invariant 
Riemann differential ideal (|2.2|) is given by the set 

R (3) {u} : = { E E ^ i3+n) fi 3+1) DlD n x u : D t 3 U = 0, 

(2.29) fl j+ V e K, n€Z+, j = W} 

and characterized by the kernel Ker D t C R^{u}, generated by the following differential relationships: 

(2.30) A/ (1) (A)-0, A/ (2) (A) = -A/«(A), 

A/ (3) (A)--A/( 2 )(A), 

where f^(X) G IC 3 {u} :— IC{u}\ D 3 u=0 , j — 0, 2, and A G K is an arbitrary parameter. The system (|2.30p 
can be equivalently rewritten in the matrix form 

/ 

(2.31) D t /(A) = <? 3 (A)/(A), g 3 (A) = -A 

\ —AO 

where /(A) :=(/(!), /( 2 ),/( 3 )(A))T G /C 3 {u} 3 . 

Now we can construct the corresponding invariant Lax differential ideal 

3 

(2.32) L 3 {u} : = {E^7 W (A) G K 3 {u} : A/ (1) (A) = 0, A/ (2) (A) = -A/«(A), 

i=i 

A/ (3) (A) = -A/( 2 )(A), 5j G /C, j = 173} , 
whose D^-invariance is looked for in the linear matrix form 

(2.33) D x f(X) = l 3 [u;X]f(X) 

for some matrix l 3 [u; X] G End M 3 . The latter satisfies the differential- functional equation 

(2.34) D t h + hD x ui = [q 3> h], 

where we put, by definition, l 3 :— l 3 [u: A], q 3 :— (73(A) G End K 3 . Making use of the derivative represen- 
tation 

(2.35) l 3 = D x a 3 , 

where a 3 := a 3 [w. X] G End R 2 , equation (|2.34l) reduces to an equivalent one in the form 

(2.36) D t a 3 = [q 3 ,a 3 ]. 
To solve matrix equation (|2.36|) . we take into account that 

(2.37) D 3 t a 3 = k 3 q 3 , Dta 3 = 

for some constant k 3 G M. Based both on (|2.37[) and on the component wise form of (|2.36[) one easily 
finds that 

Xu 2 u 3 

(2.38) a 3 [u;A] = ( Xu 2 2u 3 

— 3A 3 .x — 3A 2 ui — 2Xu 2 
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entailing the matrix 

(Xu 2iX u 3 , x 

Xu 2>x 2u 3 , x 

-3A 3 -3X 2 u hx ~2Xu 2 , x 

Thereby, having naturally extended the differential ring K, to the ring C{{a;,i}}, one can formulate the 
next proposition. 

Proposition 2.6. The Lax representation for the Riemann type hydrodynamical system 
(2.40) D t ui = u 2 , D t u 2 = u 3 , D t u 3 = 

is given by a set of linear compatible equations 





(2.41) D t f{\) = ( -A | /(A), D x f(X) 

-A 



Xu 2 . x U 3 ,: 







Xu 2 



— 3A 3 — 3A 2 mi x —2Xu 2 x 





2«3,x I /(A), 



where /(A) E Loo(K ;C ) and X € C is an arbitrary complex parameter. 

As one can easily observe, the scheme of finding the Lax representation for the Riemann type hydro- 
dynamical equations at N = 1,3 can be easily generalized for arbitrary N G Z + , that is a topic of the 
next section. 

2.4. The general case TV G Z+. Consider at arbitrary N G Z + the constraint D^u — and the 
respectively reduced invariant Riemann differential ideal (|2.2p . which is given by the set 

AT-1 

R (N) {u} : ={£ E \- {j+n) fi j+1) DiDZue1C{u}: D?u = 0. 
j=o nez+ 

J n 



(2.42) G R{{x,t}}, neZ + , j = 0,N-l}. 

The related kernel KerD t C R^{u} of the Z? t -differentiation is, owing to (|2.3j) and (|2.42|) . generated 
by the relationships 

(2.43) A/ (1) (A) = 0, A/ (2) (A) - -A/«(A), .., A/ (W) (A) = -Xf^'^X), 
which can be rewritten in a matrix form as 



/ \ 

-AO 



(2.44) 



A/(A) = <7iv(A)/(A), q N (X) = 



-A 







o '•• 
V -AO/ 

where, by definition, a vector /(A) := (/ / (2) , / W ) T G IC N {u} N := The latter 

generates the invariant Lax differential ideal 

N 

(2.45) LW{u} := {E&/ (i) ( A ) e : , 9j G /C, j = 1^, 

whose D^-invariance holds, if its linear matrix representation in the space K.( N ^{u} N 

(2.46) r>x/(A) = Zjv[u;A]/(A) 

is compatible with (|2.43[) for some matrix Zjv := In\% A] G End R N and arbitrary A £ 1. As a result we 
obtain the compatibility condition 

(2.47) D t l N + l N D x u=[q N ,l N }. 
Having represented the matrix In G End 1^ in the derivative form 

(2.48) l N := D x a N , 
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we can reduce the functional-differential equation (|2.47p to 
(2.49) D t a N = [q N ,a N } 

Now, taking into account that 
(2.50) 



D^dN = kNqN, 







for some constant G R, one can easily obtain by means of simple calculations the following solution 
to equation ()2.49j) : 







(2.51) 



a N [u; A] 
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entailing the exact matrix expression 
/ Aujv_i j2; 


l N [u; A] 
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- 1)ujv-i 
N)u N -i 
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(N - l)u N}X 



-X Nu N - 2 ,x A(l - N)u N _i tX j 



Thereby, having naturally extended the differential ring K. to the ring C{{x,i}}}, one can formulate the 
following general proposition. 

Proposition 2.7. The Lax representation for the generalized Riemann type hydrodynamical system 
(2.53) D t ui = u 2l D t u 2 = u 3l D t u N -i = u Nl D t u N = 

is given for any arbitrary N £ Z + by a set of linear compatible equations 



(2.54) D t f(X) 



D x f(X) 



where /(A) G L^p 
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/(A), 



-X z Nu N - 2 ,x A(l - N)u N - ljX J 



and X G C is an arbitrary complex parameter. Moreover, the relation- 
ships \2. 54\ ) realize a linear matrix representation of the commutator condition H.5\) in the vector space 
ICW{u} N . 



The general Lax type representation (|2.54|) . obtained above for arbitrary N E Z + , looks essentially 
simpler than those obtained before for TV = 1,4 in [31 117[ [5] and given by differential- matrix expressions 
(jl.6|) . depending at N ^ 3 on the solutions to the set of differential- functional equations (|1.7I) . As it 
was already mentioned above, it was first found in a similar form up to a scaling parameter A G C 
in [4] by means of another approach. Taking this simplicity into account, one can apply to the Lax 
type representation (|2.54j) the gradient-holonomic approach and find the corresponding bi-Hamiltonian 



THE DIFFERENTIAL-ALGEBRAIC AND BI-HAMILTONIAN INTEGRABILITY ANALYSIS OF THE RIEMANN TYPE HIERARCHY REVISITEB 



structures, responsible for the Lax type integrability of the Riemann type hierarchy of hydrodynamical 
systems ()2.53|) . 

3. The bi-Hamiltonian structures for cases N = 273 and N = 4 

To proceed with proving the related bi-Hamiltonian integrability of the generalized Riemann type 
dynamical system (|2.53j) at arbitrary N £ Z+ and with finding the naturally associated with the Lax 
representations (|2.54l) Poissonian structures, we need preliminarily to study with respect to the gradient- 
holonomic approach [21 [35J [2E] the corresponding spectral properties of the linear differential system 
(I2.46[) . Below will analyze the next cases: N = 2, 3 and N = 4. 

3.1. The case N = 2. Making use of the exact matrix expression (|2.2ip . one can define for the second 
linear equation of f|2 .28[) the corresponding monodromy matrix S(x;X) £ End C 3 ,a; 6 R, which satisfies 
[21 HE 1251 123 12H] the Novikov-Marchenko commutator equation 

(3.1) D x S(x;X) = [l 2 [u;X\,S(x;X)}, 

Since the trace-functional 7(A) := trS(x; A) is, owing to the Lax representation (|2.41l) . invariant with 
respect to the i-evolutions on the manifold M 2 , its gradient <p(x\ A) := grad7(A) 6 T*(M 2 ), x £ K, 

depending parametrically on the spectral parameter A £ C, equals 

(3.2) ip(x;\)=tT(l' 2 *S(x;\)) 
and satisfies the well known gradient relationship 

(3.3) <&(p(x; A) = z(\)r]tp(x; A) 

for some meromorphic function z : C — > C and x £ K, where i9, 77 : T*(M 2 ) —> T(M 2 ) are compatible 
[23l [25[ [27] Poissonian structures on the manifold M2, the dash sign "'" means the usual Frechet 
derivative and "* " means the corresponding conjugation with respect to the standard bilinear form on 
T*(M 2 ) x T{M 2 ). The latter naturally follows from the Lie-algebraic theory [27l [23l [25] of Lax type 
integrable dynamical systems, treating them as the corresponding Lie-Poisson flows on the adjoint space 
Q* (A) to a suitably constructed centrally extended affine Lie algebra (?(A),A £ C. By simple enough 
calculations one find that 

(3.4) V {x;\) := (^ 2 ) T = (A(S( n > - S^), S™)\ 
where we put, by definition, 

C(U) o(12 

(3.5) S(x;X) 



S (21) g(22) 

Taking into account that equation (|3.1I) can be rewritten as the system 

(3.6) = -u 2 , x S^ + 2X 2 S^ 12 \ 

^ 12) = -u 2 ] x (SW-SW) + 2\ Ul , x SW, 

= -2X 2 u 1>x (S^ -S^)-2Xu ltX S^\ 

= -2X 2 S^ + u 2 ^ 21 \ 

by means of simple enough calculations, consisting in substituting (|3 .4[) into (13.61) . one easily finds the 
gradient relationship Q3.3P in the exact differential-matrix form: 

(3.7) r? - M(M=2A(vr ^fA If" 1 

V.1 J \ (p 2 J V s u i,x u 2 , x d L +d L u 2 . x J \ cp 2 
where, by definition, z(A) := 2A, A £ C, and 



(3.8) 7?:= I _ x Q Lr? 



M _ f d- 1 u^d- 1 

d^Ul* U 2iX d~ X + 3~ l U 2lX 



constitute the corresponding compatible Poissonian structures on the functional manifold M 2 . Thus, one 
can formulate the following proposition, before stated in [7] [TBI E 13] using different approaches. 
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Proposition 3.1. The Riemann type hydrodynamical system U.ty) at N = 2 is a Lax integrable bi- 
Hamiltonian flow with respect to the compatible on the functional manifold M 2 Poissonian pair f 3.8]) . 

As a consequence of Proposition (|3.1j) easily one states that there exists (23] [27] [25] an associated 
infinite hierarchy of commuting to each other integrable bi-Hamiltonian flows on M 2 

(3.9) d( Ul ,u 2 y/dt n := -(771?- 1 u 2 , x )\ 
where t n € K, n £ Z_|_, are the corresponding evolution parameters. 

3.2. The cases N = 3 and 4. For the case N = 3 there was proved in [2 OH] that the corresponding 
Riemann type hydrodynamical system (|2.40p is a Hamiltonian dynamical system on the 27r-periodic 
functional manifold M 3 := {(ui, u%, u 3 ) J £ C^ 00 ^ (M/(27rZ); R 3 )}, whose exact Poissonian structure is 
given by the differential-matrix expression 

/ d- 1 u^d- 1 

(3.10) V=\ d^ui.x ui^ 1 + d^ 1 u 2 , x d^ 1 u 3tX 

V u^d- 1 

acting as a linear mapping ?y : T*(M 3 ) —t T(M 3 ) from the cotangent space T*(M 3 ) to the tangent space 
T(M 3 ). Making use of the exact matrix expression (|2.39[) . one can define for the linear equation (|2.33[) 
the corresponding monodromy matrix S(x; A) € End C 3 ,a; £ M, which satisfies the Novikov-Marchenko 
commutator equation 

(3.11) dS(x;X)/dx = [l 3 [u;X],S(x;X)}. 

Since the trace-functional 7(A) :— trS^rr; A) is invariant with respect to the K 9 x, i-evolutions on the 
manifold M3, its gradient f(x; A) :— grad7(A) 6 T*(M 3 ), depending parametrically on the spectral 

parameter A G C, equals 

(3.12) ip(x: A) = tr(l'*S{x;X)) = {-3X 2 S {23 \ X(S x n) + S {22) - 2S X 23 \S {21) + 2S X 32) ) T 
and satisfies the gradient relationship 

(3.13) ■dtpix; A) = z(X)r]<p(x; A) 

for some meromorphic function z : C — 5- C and all x £ R, where d : T*{M%) — > T(M%) is a second 
compatible with (|3.10l) Poisson structure on the manifold M3. In this case the standard calculations, 
consisting in substituting (|3.3j) into (|3.2p and reducing the result to the form (|3.3|) . give rise by means 
of slightly cumbersome calculations to the equivalent to (|3.13|) relationship 

(3.14) <p(x; A) = \H~ x r)(p{x\ A), 
where 



(3.15) 



du 3 + u 3 d -u 2 , x -du 2 -2u x d + 2d U3 ™ 3 1 ^ 
u 2 . x -u 2 d d Ul + Ul d -2-28^ 

1*3,3; 

-2du 1 + 2^^d 2-2^3-a 



v ~ a( — — } / 

is the second compatible with (|3.10[) co- Poissonian structure on the functional manifold M 3 . 

A completely similar results hold in the case N = 4. Namely, the corresponding Novikov-Marchenko 
equation 

(3.16) dS(x; X)/dx = [l 4 [u; X],S(x; A)], 
jointly with the expression 

(3.17) v?(z;A) = tT(l£S(x;\)) = {-4X 3 S x 2i \ -4X 2 S x 34 \ 

\(S^ + S x 22 ^ + 5i 33 ) - 3Si 44 >), S x 2r > + 2S^ + 3S^)\ 
is reduced by means of simple but cumbersome calculations to the gradient relationship 

(3.18) tp(x; A) = X 3 ^- 1 nip(x- 1 A), 
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where the Poissonian structures n and $ : T*{Mi) — > T{Mi) are, respectively, inverse to the co- 
Poissonian operators 



(3.19) 



and 
(3.20) 



V 



/ 









-d 





d 





-d 
















"4,x 


"4,1 




-3M4,£c 





«4,x 
1 "2.. 



3«4,x 



<9u 3 



-u 3 d 



113.., 

"4,1 



"4 X 



duo 



du\ + uid+ 



-a 

"4.x 

g "i,x 

"4,x 

2 rtir2j ~ 2ui,x«3,a:)5- 



2 [ U 4,x( U 2,x 



u 3 d 



u 2 d 



-du 2 - u 2 d- 

° «4,x 

du 3 



dui + u\d+ 
+^ du 2 



-u\d — u\d 



-5 + 00 



2,3! \ 



dU3 x 



5+ , 

(«1"3). 



)d- 



«4— "2"2,x Q 1 
"4.x 

qU&~ "2"2,x 



J 



on the functional manifold M4. Below we will reanalyze the bi-Hamiltonian structures of the generalized 
Riemann type hierarchy from the geometric point of view, described in detail in [251 221 13] • 



3.3. Poissonian structures: the geometric approach. To construct the Hamiltonian structures, 
related with the general dynamical system (jl.2p , by means of the geometric approach it is first necessary 
to present it in the following equivalent form: 



(3.21) 



dui/dt = u 2 — u\U\^ x 
du 2 /dt = 113 — uiu 2tX 

duj/dt = Uj + \ — U\Uj >x 

duN jdt = —U\Un,x 



> := K[u] 



where K : Mn — > T{Mji) is the corresponding vector field on the functional manifold Mjv for a fixed 
N G Z+. Then the following proposition [H [25J [28] holds. 

Proposition 3.2. Letip G T*(Mjv) be, in general, a quasi-local functional-analytic vector, which satisfies 
the condition tp ^ ip '* and solves the Lie- Lax equation 



(3.22) 



L K ip := ip t + K''*tp = grad £, 



for some smooth functional C G D(Mn). Then the dynamical system \3.21\) allows the Hamiltonian 
representation 



(3.23) 



K[u] = -77 grad H {v) [u], 



^ {n) ,K)-C 



0))' 
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where rj : T*(Mn) — > T{Mn) is the corresponding Poissonian structure, invertible on Mtf. Otherwise, 
if the operator rj^ 1 : T(Mm) — > T*(Mn) is not invertible, the following relationship 

(3.24) gia,dH M [u} = -7 ] - 1 K[u} 
holds on the manifold Mjq. 

Remark 3.3. Concerning a general solution to the Lie-Lax equation (|3.22p it is easy to observe that it 
possesses the following representation: 

(3.25) V = 4> + grad £, 
where 

(3.26) Lk^ := $ t + K' = 0, L K £ = £. 
The related co-Poissonian operator 

(3.27) rT 1 := ip'-ip 1 '* = ij>' - , 

owing to the Volterra symmetry condition (grad £)' = (grad £)''*, persists to be unchangeable. 

Based on Proposition 13.21 it is enough to search for non-symmetric functional-analytic solutions to 
(13.221) making use for this, for instance, of special computer-algebraic algorithms. 

3.3.1. The case N = 3. The corresponding dynamical system 

dui/dt = U2 — UiU± tX 1 

(3.28) duijdt = u 3 — u%U2,x / := K[u] 

du 3 /dt = -u\u 3<x J 

is defined on the functional manifold M 3 . A vector ip := (ip 1 , ip 2 > i } 3 ) 1 £ T*(M 3 ) satisfies owing to (|3.22l) 
a set of functional-differential equations 

dip l /dt + uiip l x - u 2 . x ip 2 ~ u 3.xi> 3 = S£/6ui, 

(3.29) dipz/dt + fa-iui^x = S£/5u 2 , 

dxp 3 /dt + -02 + (uiip 3 )x = S£/Su 3 

for some smooth functional £ 6 T>(M 3 ). System (|3.29l) can be slightly simplified, if to use the differenti- 
ation D t : M 3 -> T(M 3 ) : 

(3.30) D t ip 1 - u 2 , x ip 2 ~ "3,x03 = S£/5u 1 

D t tp 2 + '<Jj 1 + ui. x ip 2 = 6£/5u 2 , 
D t tp 3 + tp 2 + ui tX ip 3 = S£/5u 3 , 
which is equivalent to the scalar functional-analytic expression 

(3.31) £>t 2 Mi) = D t {a5£/8u l ) + {D t a)5£/5u 1 + (D^a)5£/Su2 + 6£/5u 3 

on the alone function ip 1 s K{u\, where we have denoted the function a := u 3x , satisfying the useful 
differential relationships 

(3.32) D t a = u\ }X a, D^a — U2, x a, D^a = 1, Df a = 0- 

The corresponding functions 2 j tp 3 G /C{u} are given by the functional-operator expressions 

(3.33) V 2 = a~ l D 1 T l (-a^ l +a5£/5u 2 ), 

ij) 3 = a~ 1 D^ 1 (— a%jj 2 + a5£/Su 3 , 

easily following from (|3.30p . In the special case £ = equation (|3.3ip reduces to 
(3-34) A 2 Mi) = 0, 

whose functional-analytic solutions can be found analytically by means of both differential-algebraic 
tools, devised in [17] . and modern computer-algebraic algorithms. In particular, making use of the 
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differential-algebraic approach of the work |17j . one can easily enough to find, amongst many others, the 
following solutions to (|3.31|) : 



2n- 



(3.35) 4> l = Ui, x /2, C( J1 ) = - I (2u 3 + u 2 u ltX )dx; 

z Jo 

-01 = ui tX u 3 - Uiit 3 , x , = 0; 



and 

(3-36) ip 1 = -u 3;X /2, £ (l5o) = 0; 



= u 2 u 3 , X: £( tfl ) = 0; 



mu2 u\ 

ipi = Ux tX u 2 - 2u 3 - ( ^-2)^3,*, £(# 2 )=0; 

W-3 3 

/ ulu 3:X 

1P1 = v>i,x — 2^"' L (^3) = °; 

Vi = > Mtf 4 ) = 0; 

u 3 



giving rise to the generating vectors 

i>( v ) = (ui, x /2, 0, -u 3 lul tX /2 + u 3 x u 2 , x ) J , £( v ) = \ J 2,r ( 2u 3 + u 2 u 1 . x )dx; 



(3.37) 



^(tf) = {Ul,xU 3 - UlU 3>x , UlU 2 ,x-Ul, x U2, 2u 2 - Hl^l.g + " 3 "i.*^" 3 " 2 ' x )T ; £ (]?) = Q; 



/ _ / ^2 ^2 a; 1 "3 "2 x U2U2 x "2"i x "3 "2 x ^3^1.x "2 x 1 

V(0) - l M 2 - "g^j- + 24^f7' 2^ 2^ft 6«L ' + 



3,x "™3,x 
3 „, .,4 „, .,6 



:!()».! 



«2«2., . U 2 U>2 ! , U 3 U 2 x »2"l,-»2. , »3«1,.«2 I . u 3«2xN T /■ _A 

1 7$ h ^ h ^r-_ 30wn) ' M 9 ) ~~ U; 



and 



(3.38) 



V'(tfo) = (~U3,x/2, u 2 , x /2, u 1 . x /2-u 3 1 x ul x /2) T , £ (t , o) =0; 

^(dj.) = C"2«3,X) -UlU3,z, UlU2M3,x% 1 - 1*^3,3;% 2 /3) T , £( #1 ) = 0; 



= («!.,«, - 2.3 - - <) U3 , X , -u 2 - ^ + agg - + u Sf , 

_ "l"l.xtt2,x 1 "i"2M3.x _ «1«2 "3, x , "3"l,x , "2"3,x _ 2tl 3"l x "2,x 

1 u 3 ,x 2u| 3m| h us,« + 15«| ^ h 

I 1 M l M 2,x _ "3"l,x"2.x _ , "3"1,.»2,, _ ^3^2^x_ \ y _ n 

+ 2« 3iI "3^7 2~^ 12u|^ 1 3Hf^ 15ug ia . ^ 'M^a) - u > 

, _ / "l"3,x "l"3,x— 6M3 (tij— 6Ujul)« 3j x . tl2-«l«l,x 1 «2«2,x(3lilll3-M \) \ T /* _ n. 

V(tf 3 ) - 2Sf~' 6S| ' 12H| 1 ^ V 3^ r >Mt> 3 ) - U ' 

^ _ / _ UlU3,i "1"2,i-"Li"2 "I"l,i-2U2 _|_ 2m 2 ,x— \ J /• Q. 

' C^4) V I) 31 li3' U3 '^3 ^3 x ' ' (^4) ' 
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As a result of (|3.23|) and expressions (|3.37[) one easily obtains the corresponding co-Poissonian operators: 



(3.39) 



d 




V 



-ui x u 3 



id 



u 3,l d 



-{u %x u 3 *d + du 2 , x u 3 l) ) 



9u 3,x 



iT 1 := M 



( u 3 d + du 3 



-2^1 + 2^- 

"3,x 



-u 2 ,x - du 2 
uid + dui 



2 - 2-^d 

"3,x 



-2 - 2d-^- 

«3,!B 

- 'IP 9 " 

tl 3«?. x +2«3«2,x ^ 



and 



(3.40) 










-du 2 , x u 3 1 



-U2,xU 3>x d 2( U ix U 3,x d + du 2,x u 3,x) 





u 2 <9 - u 2 u 3:X /u 3 

Uid - MiU3 ia; /M3 + 



uiu 3<x /u 3 - 
-u 2 u 3 , x /u 3 - dui 



+ulu 3jX /ul 
ui^du^ 1 + u 3 1 duiu 2 - 



and so on. The second expression of (|3.39p coincides exactly with the co-Poissonian operator (|3.15p . 
satisfying the gradient relationship (|3.14|) . which proves the bi-Hamiltonicity of the Riemann type 
equation (|3.28p . Concerning the next two expressions (|3.39p and (I3.40p it is easy to observe 
that they are not strictly invertible, since the translation vector field d/dx : M 3 — > T(A1 3 ) with 
components (ui )X , u 2jX , u 3 , x ) T £ T(M 3 ) belongs to their kernels, that is f] (ux jX , u 2tX , u 3 . x ) T = 
= i!)q 1 (ui. x , u 2 ,x, u 3 .x) J ■ Nonetheless, upon formal inverting the operator expression (|3.39p . we obtain 
by means of simple enough, but slightly cumbersome, direct calculations, that the Poissonian ?y-operator 
looks as 



(3.41) 



d- 1 ui^d- 1 

7] = \ d~ l Ux >x U 2 .xd^ + d~ l U 2>x d~ l U 3 , x 

u 3 . x d^ 

coinciding with expression (|3.10[) . and the corresponding Hamiltonian function equals 

(3.42) Hfr) := {i/j^,K) -£(,,)= / dx(u hx u 2 - u 3 ), 



Concerning the operator expression (|3.40p the corresponding Hamiltonian function 
(3-43) H(# o) := (ip { # o) ,K) = / dx{u 3 u 2 . x ) 



It is evident that these Hamiltonian functions (|3.42p and (|3.43p are conservation laws for the dynamical 
system (j3.28[) . which were also before found in [3]. 

Remark 3.4. It is worth to mention here that as the sum of vectors X^=o s j^($j) ^ T*(M 3 ) with arbitrary 
Sj el,j = 0,4, solves the determining equation (|3.22p . the corresponding operator 

4 

(3.44) 9 1 : ; 

will also be a co-Poissonian operator for the dynamical system (|3.28p . 
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3.3.2. The case N — 4. The Riemann type hydrodynamical equation (|1.2|) at N = 4 is equivalent to the 
nonlinear dynamical system 



(3.45) 



Ul.t 
U2.t 

U 3 ,t 



U 2 
W4 



"4,t 



-uiu itX 



:= K[ui,U2,U3,Ui\, 



where K : M4 — > T(M/x) is a suitable vector field on the smooth 27r-periodic functional manifold M4 := 
C°° (M./2irZ; R 4 ). To study its possible Hamiltonian structures, we need to find functional solutions to 
the determining Lie-Lax equation (|3.22[) : 



(3.46) Vt + K '* ip = grad C 

for if) <G T*(Mn) and some smooth functional C G D{M^), where 



(3.47) 



K 



f -du\ 


1 




— uxd 1 


-U 3 ,x 


-uxd 


\ "«4i 






\ 






1 

-uxd J 



K 



( 

1 


V 



-1*2,1 

dux 
1 
Q 



-U3,x 



dux 
1 






dux J 



are, respectively, the Frechet derivative of the mapping K : M 4 — > T(Mi) and its conjugate. The small 
parameter method [25j 13] , applied to equation (|3.46|) , gives rise to the following its exact solution: 



(3.48) 



Ux, x U 3x . 



in) 



271 



(u4Ux, x - u 2 u 3 ^ x /2)dx. 



2u 4x U 4x 

As a result, we obtain right away from (I3.23|) that dynamical system (|3. 451) is a Hamiltonian system on 
the functional manifold M4, that is 

(3.49) K = grad H M , 

where the Hamiltonian functional equals 



(3.50) 



C (v) - 



and the co-implectic operator equals 

/ 



(3.51) 



V 









-d 





d 





-d 








«4,i 




"l,x Q 
«4,l 



(uiw 4jX - u 2 u 3jX )dx 



«4,i 

M4,x 
"4,x 

2L"'4,xV"'2,z 2ux, x u 3iX )d-{ 

+d(ul x - 2ux, x U3,x)U4 x ] 



2 i u A,x( u 2,x 



J 



The latter is degenerate: the relationship ri~ 1 {ux. x , u 2>X j u 3.x, "4,a;) T = exactly on the whole manifold 
M4, but the inverse to (|3.5ip exists and can be calculated analytically. 

Proceed now to constructing other solutions to the generating equation (|3. 461) at the reduced case 
C = 0, having rewritten it in the following form: 



(3.52) 



D t 1px ~ U2,x4>2 ~ u 3.x^P 3 ~ U4.X1P4 = 0, 

Dt1p 2 +^1+ U l,xlp 2 = °! 

D t ip 3 + 4> 2 + ux, x ip 3 = 0, 

D t tf) 4 + lp 3 + Ux,xi>4 = °> 



where if) ■— (V^t "02j ^3> G T*(M 4 ). The latter is equivalent to the system of differential-functional 
relationships 

D t (atpx) ~ iftiDta ~ if) 2 D^a - — ^^Dfa = 0, 

D t (aip 2 ) + a ^x = 0, D t {aip 3 ) + aif) 2 = 0, D t (aif) 4 ) + aip 4 = 0, 



(3.53) 
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where we have denoted the function a := u 4 x . From p.53|) one ensues, by means of simple calculations, 
the following generating differential-functional equation 

(3.54) D t 2 (a^i)=0 



on the manifold M4. As above in the case N = 3, the obtained equation (|3 . 54[) can be easily enough 
solved by means of the differential-algebraic approach, devised before in [3] , which based on the following 
observation: owing to the equality D\ol = the set 

(3.55) A{u} := {f a + f x D t a + f 2 D 2 t a + f 3 D 3 t a + f 4 D*a £ JC 4 {u} : £ K 4 {u},j = C\4} 

generates a finite dimensional invariant differential ideal of the differential ring JCi{u} :— IC{u}\ D 4 u=0 . As 
a result of constructing solutions to the functional-differential equation (|3.54[) . as non-symmetric elements 
of the ideal (|3.55p , one finds the following expression: 



(3.56) 



7u 3 



u 2 u ltX + u x u 2tX + ~ (uiu 3 - ul/2) x ]y 



— 0. 



Now taking into account the expression (|3.23[) one easily obtains the second co-Poissonian operator 



^(1?) — ^(i) m ^ ne following matrix form: 



d- 1 = 



-3u4,x 



3Ui iX 



du 3 



u 3 d 



u 2 d 



-du 2 - u 2 d- 
-u 3 d^- 

° W4,x 

du 3 

dui + wi<9+ 
du 2 



"3.; 

U4 ; 



du 2 



—u\d — u\d 



-5 + d( U4 ~ U2U2 ^ ) 

«4,x B4,i 



-du 2 — u 2 d— du\ + uid+ 



+u 2 d^ 

"4, a; 

+^ du 2 



5 + (- 



U4 

(uiu 3 ) x O U 3 



)d- 



d+ 



Uj — I12"2.j 

"3,x Q _ g "3,x 
u 4,x u 4,x 



(3.57) 

satisfying jointly with (|3.51[) the gradient relationship (|3.18l) . 



4. Conclusion 

A differential-algebraic approach, elaborated in this article for revisiting the integrability analysis of 
generalized Riemann type hydrodynamical equation made it possible to construct new and more 

simpler Lax type representation for the general case N £ Z + , contrary to those constructed before in 
[17i [3] . This representation, obtained by means of the suggested differential-algebraic approach, proves 
also to coincide up to the scaling parameter A € C with that first obtained by Z. Popowicz in [4]. It is 
worth to mention that the corresponding Lax type representations for the generalized Riemann type 
hierarchy are well fitting for constructing the related compatible Poissonian structures and proving their 
corresponding bi-Hamiltonian integrability. The corresponding calculations are fulfilled in details for case 
N = 1,2 and preliminary results are obtained for cases N = 3,4 by means of the geometric method, 
devised in [25l[28j[3]. It was also demonstrated that the corresponding differential- functional relationships 
give rise to the suitable compatible Poissonian structures and present a very interesting mathematical 
problem from the differential-algebraic point of view, which is planned to be studied in other work. 
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